Probing leptoquark production at IceCube 
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We emphasize the inelasticity distribution of events detected at the IceCube neutrino telescope as 
an important tool for revealing new physics. This is possible because the unique energy resolution 
at this facility allows to separately assign the energy fractions for emergent muons and taus in 
neutrino interactions. As a particular example, we explore the possibility of probing second and 
third generation leptoquark parameter space (coupling and mass). We show that production of 
leptoquarks with masses > 250 GeV and diagonal generation couplings of 0(1) can be directly 
tested if the cosmic neutrino flux is at the Waxman-Bahcall level. 



I. INTRODUCTION 

Leptoquarks are SU (3)-colored particles which simul- 
taneously carry non-zero baryon and lepton quantum 
numbers. They are predicted in several models (such as 
SU (5) [J or Pati-Salam S77(4) 0) addressing the uni- 
fication of the lepton and quark sectors of the standard 
model (SM). In such models, the masses of the lepto- 
quarks are generically superheavy, on the order of the 
GUT scale, which puts them out of reach of direct ex- 
perimental access. Nevertheless, since leptoquarks with 
electroweak scale masses are not disallowed for any fun- 
damental reason, it is of interest to conduct experimental 
searches to delimit their properties @. It is, of course, 
important to note that in order to avoid rapid baryon de- 
cay, the simultaneous trilinear coupling of the leptoquark 
to a purely hadronic channel needs to be excluded (4J. 

In general, the couplings of the leptoquark need not be 
generation-diagonal, and the problem of extracting limits 
on couplings and masses is complicated by the presence 
of a large-dimensional parameter space. Experiments at 
HERA have placed lower limits of 0(300 GeV) on first 
generation leptoquark masses, for trilinear couplings of 
electroweak gauge strength [5|. Similar bounds, under 
the same assumptions, have been found at LEP from 
their search for anomalous 4-fermion vertices Q. For 
first generation leptoquark trilinear couplings which are 
much smaller than gauge strength (as is the case for the 
Yukawas in the SM), the mass bounds are greatly weak- 
ened. 

At the Tevatron, the leptoquarks could be produced 
in pairs, with identification made through decay topolo- 
gies. In this way, the bounds are not dependent on the 
trilinear couplings, except for decay branching fractions. 
In the case of first Q and second @] generation lepto- 
quarks, the final state topology consists of 2 hadronic 
jets + 2 charged leptons, and the resulting lower limits 
on the leptoquark mass are around 250 GeV. In the case 
of the third generation, a lower limit of 219 GeV has 
been recently reported by the D0 Collaboration by 
tagging on 2 b jets + missing energy. For this value of the 



leptoquark mass, the decay into tr is largely suppressed 
compared to the bv T channel, so that the mass bound is 
nearly independent of even the branching fraction. As 
the explored mass region becomes larger, the tr channel 
becomes more available and thus the mass limit obtained 
is pushed a bit lower (to > 213 GeV). 

In this work, we explore the possibility of probing sec- 
ond and third generation leptoquark parameter space 
(coupling and mass) with the IceCube neutrino detection 
facility This experiment, located below the surface 
of the Antarctic ice sheet at the geographic South pole, 
is required to be sensitive to the best estimates of po- 
tential cosmic ray neutrino fluxes. When completed, the 
telescope will consist of 80 kilometer-length strings, each 
instrumented with 60 10-inch photomultipliers spaced by 
17 m. The deepest module is 2.4 km below the surface. 
The strings are arranged at the apexes of equilateral tri- 
angles 125 m on a side. The instrumented (not effective!) 
detector volume is a cubic kilometer. IceTop, a surface 
array of Cerenkov detectors deployed over 1 km 2 above 
IceCube, augments the deep-ice component by provid- 
ing a tool for calibration, background rejection and air- 
shower physics. The expected energy resolution is ±0.1 
on a logio scale. Construction of the detector started in 
the Austral summer of 2004/2005 and will continue for 
6 years, possibly less. At present, data collection by the 
first 9 strings has begun. 

The event signatures are grouped as tracks, showers, or 
a combination of the two. Tracks include muons resulting 
from both cosmic ray showers and from charged current 
(CC) interaction of muon neutrinos. Tracks can also be 
produced by r leptons arising in ultra-high energy v T CC 
interactions. Showers are generated by neutrino collisions 
(y e or low energy v r CC interactions, and all neutral 
current interactions) inside or near the detector, and by 
muon bremsstrahlung radiation near the detector. 

The experimental situation is greatly simplified for 
neutrino energy E v > 10 6 GeV. A cut at this en- 
ergy is sufficient to reduce the great majority of back- 
ground from muon bremsstrahlung and tracks arising 
from muons produced in cosmic ray showers. Moreover, 
the flux of atmospheric neutrinos is low above this en- 
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ergy [ll| , so this cut generates a very pure sample of ex- 
traterrestrial neutrinos. Of particular interest here, for 
E v > 10 6 GeV, there is sufficient energy resolution (±0.2 
on a logio scale) to separately assign the energy fractions 
in the muon track and the hadronic shower, allowing the 
determination of the inelasticity distribution. Similarly, 
in the energy decade 10 6 5 < E^/GeV < 10 7 5 one can 
expect good resolution (less than 5%) in "double bang" 
events generated by incoming v^s. Again, this will allow 
a reasonable measurement of the inelasticity distribution. 

In this study, we emphasize the inelasticity (y) distri- 
bution of events as an important tool for detection of new 
physics. In particular, we will find that the y distribution 
of events generated through resonant leptoquark produc- 
tion differs substantially from the SM prediction. If the 
event rate for the new physics turns out to be comparable 
to SM expectations, then the y profile of the measured 
data can be used to probe the coupling-mass leptoquark 
parameter space. The outline of the paper is as follows: 
In Sec. [II] we derive the relevant y-distribution of events 
generated through production and decay of a scalar lepto- 
quark under the assumption of diagonal generation cou- 
pling. Armed with this distribution, in Sec. lIIII we present 
a statistical method for assessing the significance of dis- 
covery criteria. Our conclusions are collected in Sec. [TV] 



II. LEPTOQUARK PHENOMENOLOGY 

A general Lagrangian for SU(3) C ® SU(2) L <g) U(1) Y - 
invariant flavor-diagonal leptoquark couplings has been 
presented in To illustrate our proposal, we con- 

sider the simple case of SU(2)-singlet scalar leptoquarks 
Si which interact with quarks and leptons through the 
Lagrangian 

£l Q = L QiL iT2 LiL + gR l m)Si . (1) 



assume that the interaction conserves leptoquark family 
quantum numbers separately (i.e., there is no mixing be- 
tween different families) . Thus in the following subindices 
i will be dropped, and up- and down-like quarks will be 
denoted generically with U and D respectively. 

We will be considering the inclusive v P scattering 
schematically shown in Fig [T] The cross-hatched circle 
includes both resonant leptoquark production and decay, 
as well as u-channel exchange of a leptoquark leading to 
the same final state. We will not be considering couplings 
9l{b) ^ 2, since such a coupling could lead to Landau sin- 
gularities at low energies. Hence, we will assume that the 
resonant cross section largely dominates the process, and 
the narrow width 5-function approximation will be valid. 

k' 




FIG. 1: Schematic diagram of a neutrino-parton collision, in 
which a neutrino with momentum k hits a quark with mo- 
mentum p — xP giving rise to a secondary charged lepton 
and quark with momenta k! and p' , respectively. Here x is 
the fractional energy of the struck parton in the nucleon hav- 
ing momentum P. 



Let us assume that an incoming neutrino collides with 
a proton target with center-of-mass energy s. If the neu- 
trino hits a down-like parton D, the inclusive cross sec- 
tion for the process shown in Fig. [T] in the parton model 
is given by 



Here Qi — (u,i di) T and Li = (z^ li) T stand for quark and 
lepton SU(2) left-handed doublets, um and Ur are right- 
handed singlets, and and gR are the corresponding 
coupling constants. Subindices i, running from 1 to 3, 
label the quark or fermion family. For simplicity we will 



da LQ {vP -> l-X) 
d 3 k'/E' 



dx D(x) 



daLQivD -> l-U) 
d 3 k'/E' 



(2) 

where D{x) is the corresponding parton distribution 
function (pdf). The neutrino-parton cross section reads 
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da LQ {vD -* l~U) 
d 3 k'/E' 
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(3) 
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where T = 2xs = 2s is the invariant flux, my is the 
mass of the outgoing up-like quark, and E' is the lab 
energy of the outgoing charged lepton. In the resonant 
approximation, the amplitude for the production of a left- 



handed charged lepton is given by 



l L (k')Ul(p')Dl(p) VL (k) 



(4) 



s - M 2 - iTM 
where M and T are the mass and width of the leptoquark 
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A similar expression, replacing g 2 L — > 9l 9r, holds for the 
decay through the right-handed channel. It is easy to see 
that there is no interference between DJIl and DJIr. 
Now in the narrow resonance approximation one has 



1 



(s - M 2 ) 2 + (TM) 2 



TM 



5(s-M 2 ) 



(5) 



Then, after summing over spins of outgoing fermions one 
arrives to 



1Mb! 2 ) 



91 (9l + 9r) 



s (s — mfj) 
MT 



5(s - M 2 ) ,(6) 



where we have neglected both the mass of the down-like 
quark and the mass of the outgoing charged lepton. We 
have kept instead the mjj dependence, since it is relevant 
in the case of the third family, where the outgoing quark 
would be a top. As a further assumption, we will consider 
that the leptoquark width is dominated by the Ul and 
Dv quark-lepton channels. This leads to 

r = ^{^[(i-^f+^+^a-At/) 2 }, (?) 



where Xu = mfj/M 2 
into Eq. (2) yields 



Substituting Eqs. ©, © and © 



da 



(a) 
LQ 



9l(9l+9 2 R )a-^u) 



dx 



D(x) 



<Pk>/E> " ' (gl + g 2 R )(l - Xu) 2 + g\ J ' 2xs 

x M 2 5{Q 2 - 2m p y Ey x + mjj) 5(xs — M 2 ) 
J_ 9l(9l+9 2 R )(l-^u) 

(8) 



j 2 L + g 2 R Kl-Xu) 2 + 9l 
x D(M 2 /s) 5{Q 2 -yM 2 + mfj) , 



where Q 2 — — (k — k') 2 , and the inelasticity y is defined 
as y — (E v — E')/E v , E v being the lab energy of the 
incoming neutrino. Indices a — 1, 2, 3 correspond to the 
up-like quarks U = u,c, t, respectively. After adequate 
changes of variables and integrations [l3| we find 



da 



(«) 
LQ 



dy 



* gligl+gPd- Xu) D(M 2 /s) 
2 (ffi+^Xl-Ac/) 2 



9l 



(9) 



The inelasticity y lies in the range Xu < y < 1. Note 
that the y distribution of the resonant process is approxi- 
mately flat (at the energies of interest, the Q 2 dependence 
of the pdf can be neglected), in contrast to that charac- 
terizing SM charged current (CC) processes in which 



aa SM 
dy 



2 G| TYlpEy 



Q 2 + M 2 V 



x / dx[xq(x, Q 2 ) + xq(x, Q 2 )(l - y) 2 } , (10) 



where Gf = 1.16632 x 10~ 5 GeV -2 is the Fermi constant, 
Mw is the mass of the W gauge boson, and q(x) and q(x) 
stand for combinations of quark and anti-quark proton 
pdf's, respectively The y dependence of the SM 

cross section is shown in Fig. 6 of Ref. [l5[ . In the next 
section we exploit the differing y dependences of the lep- 
toquark and SM interactions to constrain the parameter 
space of the new physics. 



III. SENSITIVITY REACH AT ICECUBE 

To evaluate the prospects for probing leptoquark pro- 
duction at IceCube, one has to estimate the "beam lu- 
minosity" , i.e. the magnitude of the (yet to be detected) 
neutrino flux. We know that cosmic accelerators produce 
particles with ener gies in excess of 10 11 GeV (we do not 
know where or how [161]), and a neutrino beam is expected 
to come in association with these cosmic rays [17j | . How- 
ever, given our ignorance of the opacity of the sources, 
it is difficult to calculate the magnitude of the neutrino 
flux. The usual benchmark here is the so-called Waxman- 
Bahcall (WB) flux 



E, 



"WB 



6 x 1CT 8 GeV cm" 2 s" 1 sr" 1 (11) 



(all flavours), which is derived assuming that neutrinos 
come from transparent cosmic ray sources [18] , and that 
there is adequate transfer of energy to pions following pp 
collisions. Here we will rely on this expression to estimate 
the event rates needed to quantify the IceCube sensitivity 
to leptoquark production. However, one should keep in 
mind that if there are in fact "hidden" sources which 
are opaque to ultra-high energy cosmic rays, then the 
expected neutrino flux will be higher [l9[. Moreover, if 
the extragalactic cosmic rays begin to dominate over the 
Galactic component at energies as low as ~ 10 9 GeV, as 
suggested recently [2(j, then the required power of the 
extragalactic sources will increase by a factor of ~ 2, 
implying a concommitantly larger neutrino flux (2lj . 

IceCube is sensitive to both downward and upward 
coming cosmic neutrinos. However, to remain conser- 
vative with our statistical sample, here we select only 
downward going events. To a good approximation, the 
expected number of such events at IceCube is given by 



M = 2vrn T T 



J dE v a tot (E u ) c^ B (E u ) , (12) 



where tit is the number of target nucleons in the effec- 
tive volume, T is the running time, and (Jtot{E v ) is the 
total neutrino-nucleon cross section. In our analysis we 
are interested only in CC contained events, for which an 
accurate measurement of the inelasticty can be obtained. 
The IceCube's effective volume for (background rejected) 
contained events is roughly 1 km 3 [22[ , which corresponds 
to n T ~ 6 x 10 38 . 

Hereafter we focus on neutrino energies in the range 
10 7 < Eu/GeV < 10 7 5 , where the background from at- 
mospheric neutrinos is negligible, but the extraterrestrial 
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' 200 220 240 260 280 300 320 340 

M (GeV) 

FIG. 2: Sensitivity reach of IceCube (90% CL) to probe 
second generation SU(2)-singlet scalar leptoquark parameter 
space (coupling and mass). For comparison, the existing limit 
(95% CL) reported by DCCollaboration [j| is also shown. 



flux is expected to be significant. Thus, we will consider 
a medium energy (E„) = 10 7 25 GeV. At production, the 



WB flux has flavor ratios v a 



2:1:0, but this 



quickly transforms to 1 : 
tions [23|. One has then 



1 : 1 through neutrino oscilla- 



6 x 1CT 23 GeY -1 cm" 2 s" 1 sr" 1 



(13) 



for each flavor a. Now it is possible to increase the ratio 
signal/background events by performing a cut in the in- 
elasticity y. Given the dependence on y of the Cg^f cross 
section, the flat behavior of the <jlq cross section, and 
the available phase space for quark production, it is con- 
venient to consider events with relatively large values of 
y. We choose here events in the range y > 0.5. With this 
cut, the integration of Eq. (fTo)) leads to 



y>0.5 



8 X 10 



" 34 cm 2 



(14) 



One has to take into account that this result carries a 
systematic error of about 20% [TJ, [23 due to uncertain- 
ties in the extrapolation of the pdf 's [25| • For illustrative 
comparison the second generation leptoquark cross sec- 
tion in the high y region is calculated from Eq. ([9]) to 
be 



vlq{(E v ))\ 



y>0.5 



2 X 10 



" 33 cm 2 



(15) 



for fiducial values ql = 9r = 1 and leptoquark mass 
M = 300 GeV. From Eqs. (T3|) and (T3J) we can now 
easily estimate the number of expected SM background 
events during the lifetime of the experiment. Taking T = 



15 yr, for each neutrino flavor a one has 

A/B (a) * 2nn T T ag°((E v ))\ y > 05 ^{(E v )) AE V 
= 2 , (16) 

where AE V = 10 75 GeV - 10 7 GeV ~ 2.2 x 10 7 GeV. 
(Note that the background is cosmic neutrinos! Todays' 
signal, tomorrow's background.) In the same way, the 
number of signal events As will be approximately given 
by 



A" s (a) ^2« T T <TL Q ({E v ))\ y > 0B ^ B ((E V )) AE V , 

(17) 

which for the above proposed leptoquark interaction is 
just a function of the couplings <?_r.l and the leptoquark 
mass M. 
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FIG. 3: Sensitivity reach of IceCube (90% CL) to probe third 
generation SU(2)-singlet scalar leptoquark parameter space 
(coupling and mass). For comparison, the existing limit (95% 
CL) reported by D0Collaboration is also shown. 

It should be noted that the neutrino induced events 
do not constitute the sole background. As mentioned 
in the introduction there are muons (produce in the at- 
mosphere) which traverse the detector and may deposit 
energy through bremsstrahlung radiation. In our energy 
bin, one may expect 10 muon traversals in 15 yr. How- 
ever, our inelasticity cut will completely eliminate this 
source of background, because of the negligible probabil- 
ity for muons to radiate 50% of their energy. 

To determine the bounds for leptoquark production, 
let us assume that 2 z/ a -events are in fact observed with 
y > 0.5. Then, at 90% CL, there will be an upper bound 
on signal events given by Ag"" 1 < 3.91 [26|. For simplic- 
ity, we consider the left-right symmetric case in which 
9l = 9r- Then, after numerical evaluation of the lep- 
toquark cross sections, the upper bounds on A/"g can 
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be translated into contours of constant likelihood in the 
M-g L plane. Our results are displayed in Figs. and 
where we show the sensitivity reach of IceCube together 
with the existing limits from D0Collaboration [1, @. 

In the case of the third family, it can be seen that the 
sensitivity is maximal for leptoquarks of M ~ 245 GeV. 
For lower leptoquark masses (in the narrow resonance 
limit) the allowed inelasticity range — and thus the lep- 
toquark cross section — becomes reduced due to phase 
space suppression, owing to the large mass of the top 
quark. 

In order to estimate the significance of the assumption 
9l = 9r, we have also considered the case of purely left- 
handed leptoquark currents, i.e. <?_r = 0. By looking 
at Eq. ((9J, it can be seen that this implies an average 
reduction in the leptoquark cross section by a factor of 
about 0.75 and 0.65 for the second and third families, 
respectively. 

IV. CONCLUSIONS 

In this paper we have introduced the measurement of 
inelasticity as a powerful tool for probing new physics in 
cosmic neutrino interactions. As an illustrative example, 
we have discussed the possibility of detecting leptoquark 
production at the IceCube neutrino telescope [27]]. We 
estimated the expected event rate at IceCube to be com- 
parable to the one predicted for cosmic ray facilities that 
make use of the atmosphere as the detector calorime- 
ter [28| . However, the ability of IceCube to accurately 



measure the inelasticity distribution of events provides 
a unique method for SM background rejection, allow- 
ing powerful discrimination of resonant processes: we 
have shown that production of leptoquarks with masses 
> 250 GeV and diagonal generation couplings of 0(1) 
can be directly tested at the Antarctic ice-cap. 

In closing, some comments are in order. First, we have 
not taken account of any systematic considerations con- 
cerning the detector - these are beyond the scope of the 
present work. Second, for reasons of simplicity, we have 
not included upcoming events close to the horizon. Al- 
though these Earth-skimming neutrinos have the poten- 
tial of nearly doubling our signal event background (and 
thus nearly halving the required observation time scale) , 
their proper consideration will require a full Monte Carlo 
simulation. 
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